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We study light propagation and gravitational leasing in scalar-tensor theories of gravity 
by using a static, axisymmetric exterior solution. The solution, which has been studied 
by Tsuchida and Watanabe, has asymptotic flatness properties and moreover is reduced to 
Voorhees’s one of Weyl’s series of prolate solutions in the case of a constant scalar field. We 
first obtain an asymptotic form of the solution near the spatial infinity in order to clarify 
the physical significance of three model parameters found in the solution. It is shown that 
an amplitude of the scalar held, a non-sphericality of the spacetime symmetry and a mass- 
like parameter in the Einstein frame are represented in terms of these three parameters. 
As for the spacetime structure, we also give brief discussion on the directional-dependent 
properties of the spacetime singularity. We then study null geodesic equations and Sachs’s 
optical scalar equations on the equatorial plane in order to investigate deflection and shear 
of light rays. Our studies are done by using a technique of the conformal transformation 
such that their results are independent of details of scalar-tensor theories of gravity owing to 
the conformally invariant properties of null geodesics. For some specific values of the model 
parameters, we analytically obtain a deflection angle of the light path and find that it can 
become negative. The appearance of a negative deflection angle indicates “reflection” of a 
light path, and we investigate under which conditions the light reflection occurs. As for the 
optical scalars, we And that the Weyl source-term shows significantly different properties 
when it is compared with that in the Schwarzschild spacetime. We therefore classify a 
space of the model parameters into four distinct regions on the basis of the qualitative 
properties of the Weyl source-term and find a close relationship between this classification 
of the parameter space and the occurrence of the light reflection. We Anally solve the null 
geodesic equations and the optical scalar equations numerically. We find that a picture of 
the thin lens is applicable and give a simple analytic model for the optical scalars. As for 
the properties of gravitational lensing, the deflection angle and the image distortion rate are 
obtained as functions of the impact parameter. Again, we find a close relationship between 
their qualitative properties and the classification above. 


§1. Introduction 


Scalar-tensor theories of gravity have been studied by many theoretical physicists 
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as natural alternatives to general relativity since the pioneering work by Brans and 
DickeP In particular, a theory of dilaton gravity has been of interest as an effective 
theory of the superstring theories at low energy scales!^ In these theories, the gravity 
is mediated not only by a tensor field but also by a scalar field. Several theoretical 
predictions in the scalar-tensor theories have been obtained,and it has been 
found that a wide class of the scalar-tensor theories can pass all the experimental tests 
in the cases of weak gravitational fields. In the cases of strong gravitational fields, 
however, it has been found that the scalar-tensor theories show different aspects 
of the gravity in contrast to general relativity. For example, it has been shown 
numerically that nonperturbative effects in the scalar-tensor theories increase the 
maximum mass of an isolated system such as a neutron star.^)“^^ In these works, 
numerical methods play an important role, and a spherical exterior solution, which is 
a static, spherically symmetric exterior solution outside a spherical body, is matched 
to the numerical solution in interpreting its results. A usual astronomical object, 
however, has a non-spherical form and is in rotation, and not only a spherical exterior 
solution but also axisymmetric exterior solutions may therefore play an important 
role in discussing the crucial difference between general relativity and scalar-tensor 
theories of gravity. 

Among the axisymmetric exterior solutions, ones corresponding to the series of 
solutions of Tomimatsu and Sato in general relativit}^^ must be of significant interest. 
It is, however, difficult to obtain their explicit forms, and the case is similar even for 
solutions corresponding to Kerr’s solution in general relativity. Instead of studying 
such the Kerr-like solutions, Tsuchida and WatanabJ^ have investigated the Weyl- 
like solutions, namely the solutions which are reduced to those of Weyl’s series in the 
case of a constant scalar field. Among these Weyl-like solutions, we have a special 
solution, which we shall refer to as the scalar-tensor-Weyl solution. This solution is 
very notable for the reason that it has asymptotic flatness properties and moreover 
is reduced to Voorhees’s solution,!^ namely a solution of Weyl’s series of prolate 
solutions, in the case of a constant scalar field. Several geometrical properties of the 
scalar-tensor-Weyl solution, especially the properties of light propagation, have been 
investigated by analytical approaches in Ref.7. On the basis of this work, we will 
study light propagation and gravitational leasing on the scalar-tensor-Weyl solution 
by numerical approaches because the local properties of null geodesics must tell us 
important geometrical information of the spacetime. 

In this paper, we first study the global structure of the scalar-tensor-Weyl solu¬ 
tion. We obtain an asymptotic form of the solution near the spatial infinity in order 
to clarify the physical significance of three model parameters found in the solution. 
It may be interesting that the Schwarzschild-like coordinates can be naturally in¬ 
troduced in doing that. It will be shown that these parameters are related to an 
amplitude of the scalar field, a non-sphericality of the spacetime symmetry and a 
mass-like parameter in the Einstein frame. The directional-dependent properties 
of the spacetime singularity will be also shown. After giving the several analytical 
results as briefly summarized in Abstract, we numerically solve null geodesic equa¬ 
tions and Sachs’s optical scalar equations.^Sl By using a technique of the conformal 
transformation, we then study deflection and shear of light rays on the equatorial 
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plane in the manner snch that the obtained results are independent of details of 
scalar-tensor theories of gravity. We numerically obtain the deflection angle and 
the image distortion rate as functions of the impact parameter and find a close re¬ 
lationship between their qualitative properties and the classification of the model 
parameter space according to the Weyl source-term. Several analytic results shown 
in the previous papers are summarized in Appendices.^ 

Throughout this paper, we use the system of units such that c = G = 1. 


§2. The scalar-tensor-theory of gravity 


We consider the simplest scalar-tensor theory in which gravitational interactions 
are mediated by a tensor field, g^i,, and a scalar field, Hereafter, a hat ( ~ ) 

is used to denote quantities and derivatives associated with 5 ^,^. An action of the 
theory is given as 


S 



uj{^) 




y/^d'^x + 

^matter 



( 2 - 1 ) 


where u;{^) is a dimensionless function of (j) and specifies the theory, represents 
matter helds, and Scatter is an action of the matter fields. The scalar field, cj), plays 
the role of an effective gravitational constant, G, as G ~ !/(/>. Varying the action by 
the tensor field, and the scalar field, (j), yields, respectively, the following field 
equations, 


r - —T J- 

^fj.iy — ^ ~r 




2 y , (p,a(p,p 


1 


+-f 


V" 


3 -|- 2Lo{(j)) 


SttT- 


duj{(p) 

d^ 


.a(P,P 


( 2 - 2 ) 

(2-3) 


where and denote, respectively, a covariant derivative and an energy-momentum 
tensor associated with Now we perform the following conformal transformation 
to a new frame called the Einstein frame with the metric, g^m, defined by 

9n>^ — i}p)9n^i ( 2 ' 4 ) 


such that 


A\g>) = 


1 


1 


/ din A((^) \ _ 


\ dip ) 


= a {ip), 


(2-5) 


4> 3 + 2oj{4>) 

where A{ip) and Oi{ip) are referred to as a coupling function and a coupling strength, 
respectively. Then the action is rewritten as 


S 


1 

Ibvr 


{R - 2g>^‘'ip^^ip^i,)^/^d'^ 

^ *^matter W.n,A\ip)g^ 




(2-6) 
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Varying the action by and ip yields, respectively, 

Gfxi/ = + 2 — —§111/9 ^‘P,aP’,/3^ ) {‘^G) 

= -4:Tra{ip)T, (2-8) 

where Vq and T^u denote, respectively, a covariant derivative and an energy-momentum 
tensor associated with g^i,. A relationship between and is given by 




2 6S^^u.,['P^,A'^{ip)g^ 


flU\ 


The conservation law for is given by 


= A^{ip)f>^^. 


(2-9) 


= a{ip)TV>^ip. 


( 2 - 10 ) 


§3. Geometrical properties of the scalar-tensor-Weyl solution 


We investigate the geometrical properties of the scalar-tensor-Weyl solution.^ 
The reduced field equations and the derivation of the solution are given in Appendices 
A and B. 


3.1. The scalar-tensor-Weyl solution 

The scalar-tensor-Weyl solution is one of the static, axisymmetric vacuum solu¬ 
tions in the Einstein frame with the metric, 




X -|- 1 


2 1 \ ^ 
X — 1 


X -I- 1 


X 


x2 _ y2 


x^ — 1 \ J 


{x -y) 


and the scalar field, 


d ( X — 1 


(32) 


where a positive constant, cr, is a unit of length, and 5, (po and d are integration 
constants. A parameter. A, is defined by 


A^ = 6‘^ + d‘^, A>0. 


(3-3) 


Now we show that the above solution is reduced to the previously known ones 
when we take specific values of the parameters, 5 and A. When A = 6, i.e., d = 0, the 
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metric is reduced to Voorhees’s prolate solutioiJ^ which is one of Weyl’s series of the 
solutions in general relativity. Voorhees’s prolate solution contains the Schwarzschild 
solution as a specific case, <5 = 1, and the Schwarzschild coordinates are related to 
the coordinates, x and y, as 


X = -1, y = cos 9, m = a, 

a 


(3-4) 


where m plays the role of a usual mass parameter of the Schwarzschild spacetime. 
When Z\ = 1, we introduce the Just coordinates by 


X = — — 1, y = cos 6, m = 6a. (3'5) 

a 

Then (I3T|) is reduced to a spherical exterior solution. 



This result tells us that the parameter. A, has the physical significance of a non¬ 
sphericality of the spacetime symmetry. The Schwarzschild coordinate, r, is related 
to the Just coordinate, y, as 


i-a 

/ 2m\ 2 

When A = 5 = 1, the spherical exterior solution (ESI) is reduced to Schwarzschild’s 
one, and, in this case, the scalar field becomes a constant. A relationship among 
these solutions is shown in Figure 1. 


A—- 5 

Scalar-Tensor-Weyl (0 , A , 5 ) - 


Weyl (0,5) 


A — 1 


5 ^ 1 


Spherical Exterior Solution (0,5 ) 


5 


1 


Schwarzschild ( 0 = m ) 


Fig. 1. A schematic depiction of series of the scalar-tensor-Weyl solutions. 


3.2. An asymptotic form of the scalar-tensor-Weyl solution near the spatial infinity 
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Now, in order to clarify the physical significance of the parameters, A and 6, we 
study an asymptotic form of the scalar-tensor-Weyl solution near the spatial infinity. 
We first define the asymptotically Schwarzschild-like coordinates, {r,9), by 



+ v‘^ + 


1 - Z \2 


+ 



1-Z\2 


2 


4^2 ^2 


(3-8) 


Z\ 


+ v'^ + 


I- A"^ 


+ 


1 - 


Av?v‘^ 


u = — -X — 1 , V = cos 6, m = 5a, 
mA 

m m / f f 

X 5 J 2r^ 


(3-9) 

(3-10) 

(3-11) 


which are reduced to (EiSI) when A = 1. The equations in the above complicated 
forms are obtained under the following requirements. 


(a) The equations dSIHl) and (Eini) are the same as those in Voorhees’s prolate 
solution, and one has X = f for A = 5 such that the asymptotic form of the 
solution becomes that of Voorhees’s one.l^ In particular, an asymptotic form of 
gtt for Z\ = (5 becomes 


gtt = -1 + 


2m 


2 

3 


, 52-1 

52 




(3-12) 


(b) When Z\ = 1, the asymptotic form of the solution is the same as that of the 
spherical exterior solution, and an asymptotic form of gu and the relation between 
the Schwarzschild coordinate, r, and the Just coordinate, Xj are, respectively, given 
by 


2m 1 

gtt = - 1 - X 

r 3 


m 


m 


= -+ 1 - 


X 


1 


1 - 


m 


J 2 


m 


+ 0 


m 


+ 3- 


1 \ m^ 


5 2r3 


+ 0 


m 

TT 


(3-13) 

(3-14) 


(c) For any values of A and 5, asymptotic forms of gu and gee become, respectively, 

, 2 \ /,^3' 


gee = r 


/m\ 

1 + 0 - (— +0 


\ r J 


m‘ 


2m 

= -IH-^0 • 

r 


m 


+ 0 


+ 0 




(3-15) 

(3-16) 


With (l3-8jl ^ (l3Tll) . asymptotic forms of the metric and the scalar field near the 
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spatial infinity are, respectively, given by 
ds‘^ = — 


2m f 1 ( \ 2 ( - 1 


3 2„ 1 

- cos o'- 

2 2 


+ 


2 m ( Z\^ 


m 


-1 


dr‘^ _|_ ^2 ^^^2 ^ 


and 

d m m? 

~i 9^o) ^ 


Z\2 _ 1 
3(52 


- COS - 

2 2 




r»3 


m 


dt^ 


(3-17) 


(3-18) 


Since {d/6)‘^ = A^ — 1, these asymptotic forms include the parameters, A and (5, 
in the forms, A^ jd’^ — 1 and (Z\^ — l)/(5^, which can be interpreted as a squared effec¬ 
tive amplitude of the scalar field and a non-sphericality of the spacetime symmetry, 
respectively. It is noted that the non-sphericality appearing in gti can be read as a 
quadrupole moment of the gravitational potential in Newton’s theory. 


3.3. Directional properties of the spacetime singularity 

The scalar-tensor-Weyl solution has a coordinate singularity at x = 1, and it 
is generally a naked spacetime singularity whose qualitative features are similar to 
those of Voorhees’s prolate solution. A spatial topology of the singularity at x = 1 
is a sphere for Z\2 ^ 1 and becomes a point for Z\2 = 1 . When Z\2 = <5 = the 
coordinate singularity at x = 1 becomes an event horizon, r = 2m, and a true 
spacetime singularity at r = 0 is hidden by the event horizon. When Z\2 ^ 1 , two 
poles, namely y = 1 and y = — 1 , on the singular surface at x = 1 have directional 
properties as shown below. 

In the scalar-tensor-Weyl solution, a scalar polynomial of the curvature tensor 
can be calculated as 

= ^ (^ - {X + l)-2(-5+/i^+l) (^2 _ y2^2zi2-3 ^ 

[-4Z\^(52 (1 - y 2 ) + 4^4^^ _ y2^^ _ 3^4 (^^2 _ ^ 2 ^ 

-|-8Z\^(5^x (l — — 2A^d‘^ (—4 -|- 3y^ — 6y^x^ -|- 7x^) 

—4A^dx (l -h y^ — 2 x^) — 7d^ (x^ — y^) 

—8d^x [l + y‘^ — 2x^) — 4d^ (l — y^ — 3x^ + 3x^)] . (3T9) 

As a simple example, we consider a path on the (x,y)—plane given by y = x — 
(x — 1)“, where o: is a constant such that 0 < (a < I. This condition on the index, 
a, ensures that inequalities, x > 1 and y < 1, are satisfied along the path. Then one 
finds that the polynomial (lOTl) along the path is expressed as the product of some 
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regular function of (x, y) and the following factor, 

{X -= (. - D'g 
K = (2A‘ - 3) q - 2 {4^ - (5 + 1) . 

While one finds that the polynomial diverges along the path when K < 0, 

one does not meet any spacetime singularity along the path when K > 0. One 
notes that the inequality, — 5 + 1 > always holds because of the definition of 
Z\: A = ^6“^ + d? > |(5|. The sign of the index, K, therefore depends on the sign of 
the coefficient, 2Z\^ — 3, as follows; 

(a) The case of 2Z\^ — 3 < 0 

The index, K, always becomes negative for any positive value of a, accordingly 
the pole, {x,y) = (1,1), is a spacetime singularity. 

(b) The case of 2Z\^ — 3 > 0 

When a > ao = 1 — {25 — 5) / (2Z\^ — 3) > 0, ill > 0 holds, and one does not meet 
a spacetime singularity at the pole, {x,y) = (1,1), along the path. 

We now summarize the conditions on the parameters, A and 5, such that the 
index, K, can be positive. When the inequality, 5 < 2/5, holds, the parameter, uq, 
defined for the above case (b) always becomes greater than unity. Since 0 < a < 1, 
two conditions, A'^ > 3/2 and <5 > 5/2, must hold in order for K to be positive. 
From the another point of view, one can think of ao as a function of A and 5. 
Under the above conditions, > 3/2 and 5 > 5/2, ao is a monotonic function of 
A and 5 and has a minimal value as oq = (l + \/l9) /6 for A = 5 = (5 + \/l9) /2. 
Consequently, the index, K, is always negative for any values of A and 5 in the case 
of 0 < a < (1 + v/l9) /6. 

§4. Analytic results on null geodesics and optical scalars 


We study the properties of irrotational null geodesics and optical scalars in 
the scalar-tensor-Weyl solution by analytic approaches. Hereafter, we consider the 
parameter, 5, to be non-negative such that m = 5a > 0 holds. Accordingly, we have 
0 < 5 < A in scalar-tensor theories of gravity and 0 < <5 = A in general relativity. 

4.1. Null geodesics 

Because of the conformally invariant nature of null geodesics, we consider null 
geodesics on the spacetime with the metric, g^y = A~‘^g^y. Hereafter, a hat denotes 
geometrical quantities associated with g^y, and the corresponding geometrical quan¬ 
tities associated with g^y are denoted without a hat. Let = dx^/dv be a null 
geodesic tangent associated with g^y^ where v is the affine parameter. Then a null 
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geodesic tangent, = dx^/dX, associated with is given by 

rj yv, 

dX = A-^dv, F = — = A^k^^. 

dX 

The null vectors, k^ and k^, satisfy the following geodesic equations, 

= 0 . 


(4-1) 


(4-2) 


Since it is difficult to study generic null geodesics analytically, we consider spe¬ 
cific null geodesics which are always on the equatorial plane, y = 0. In this case, the 
geodesic equations are reduced to 


W^ = l- 


h^ 


X — 1 


t = 


cr^ — 1) + 1 

X + 1 \ : h 


X — 1 


2(5 


= 1-V{x), 

X — ^ 


a‘^{x‘^ — 1) yx + 1 


(4-3) 


where a dot denotes a differentiation with respect to the affine parameter. A, a 
constant, h, is an impact parameter and IT is a function of x defined by 


dx 1 
dW^a 



(4-4) 


One of the particular orbits of interest is a circular orbit and it is determined by 
the conditions that IT = 0 and dV/dx = 0. These conditions are reduced to x = 25 
when 6 > 1/2, and the circular orbit cannot exist when 6 < 1/2. As will be shown 
later, 5 = 1/2 is a critical value such that many geometrical quantities for 5 < 1/2 
show very different properties compared with those for 6 > 1/2. 

In the asymptotic region such that r » m (see Section 3.2), the geodesic 
equation can be approximated by 


dr 


2 






(4-5) 


which is compared with the following exact result in the Schwarzschild spacetime. 


f .^V - 1 - — 

\dx) 



(4-6) 


With the reduced null geodesics equation (ITHl) . a deflection angle, a, of the light 
path is obtained as 



x^ - 1 




' x^ — 1 / X -|- 1 


X2-1 


X — 1 


2(5 


Xq + I 
Xo - 1 


-2(5 


:dx — TT, 


- 1 


(4-7) 
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where a parameter, Xq, is related to the impact parameter, h, by 



(4-8) 


Since a deflection angle can be defined only for the scattering orbit, the allowed 
range of Xo is given by Xo > 25 for 5 > 1/2 and Xq > 1 for 5 < 1/2, accordingly the 
possible range of the impact parameter is given by 


/i > /i* 


2m\l 1 - 

Am 

0 


25 + 1 
25-1 


(5 > 1/2) 

(5 = 1/2) 
(5 < 1/2), 


m = 5 it > 0. 


(4-9) 


Specific values of /i*/2m are, for example, 2.161, 2.543, 2.598, 2.667 and 2.690 for 
5 =0.53, 0.7, 1.0, 1.5 and 2.0, respectively. The coordinate, x, decreases to the 
“perihelion”, x = Xq, and then increases along the scattering orbit. 

We obtain an analytic form of the deflection angle, a{xo), in the following specific 
cases. 


(a) 5 = 0, Z\2 = 2iV (IV = 1, 2,3, • • •) 


a(xo) = 




I r r ^ i-in^p - i)!!(iv -1)! \ 

I V Xo'^ — 2P{N — p — l)\{pl)‘^Xo^P J 

/ IV-1 3(iV-l)(lV-2) 

V ~ 2 xo^ 16x7 

-TT (Xo ~ 1) 

7rZ\^ 



(4-10) 

One notes that it is the case of m = 0. The deflection angle is a negative, 
monotonously increasing function of Xq for any value of N, as is depicted in Figure 
2a. In addition, the deflection angle can be regarded as a monotonously decreasing 
function of N. 


(b) 5 = 1/2, = 4L + 2 (L = 0,1,2, • • •) 

(i) L = 0 


a{xo) 




TT 


Xq + I 

y^Xo{Xo + 2) 


— sin 

TT 


1 

Xo + 1 


873 




> 0 

1 


{Xo ~ 1 ) 
(Xo >> 1) 


+ 1 



(4-11) 
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(ii) L > 1 


a{xo) 

Il{xo) 


Xq + I 
Xq {xo 2) 
(2L)! 


4L+1 


■^piXo) 


Xo{Xo + 2) 

L 2n 

EE (2L — 2n)!(2n — k)\k\ 

LX —0 k —0 

, _ 

^ ^ (2L - 2n - l)!(2n -k + l)\kl 

(2p- 1)!! /vr 1 


Il{xo) - vr, 

V{Xof^^~"'^~^CiXo)^Be{Xo), 


i2p)\\ 


+ sin 

2 Xo + 1 


\/Xo{xo + 2) \ (2p — 2r — 2)!! 


1 


Bpixo) 


+ 1 ^ {‘2‘P — 2r — 1)!! (Xo + l)2p-2r-l j ’ 

{2p)\\ ^Xo{xo + 2)^{2p-2r-l)\\ 1 


(2p+l)!! Xo + l ^ {2p — 2r)l\ (xo + 1)^P“^^’ 


where ??(xo) = 2/(xo + 1), ^(a:o) = 1 — 2xo — Xq^ and i = L + n — k. 
(c) (5 = 1/2, Z\2 = 4L + 4 (L = 0,1, 2, • • •) 

-i 4L+3 


(4-12) 


a 


(xo) = -2v^Xo(xo + 2) 


L 2n+l 

Jl(Xo) = X] 


Xq+I 

Xq{xo 2) 

(2L + 1)! 


n=0 /c=0 
L 2n 

+EE 

n=0 fc=0 


{2L — 2n)\{2n — k + 1)!A:! 


(2L + 1)! 


{2L — 2n + l)!(2n — k)lkl 


Jl{xo) - TT, 

V{Xof^^~'^'^^{Xo)’'Ae+l{Xo) 

p{xof^^~'^'>'^^C{xo)’'Be{xo), 


(4-13) 

where £ = L + n — k. The functions, r]{xo), ^(xo), Ap(xo) and Bp{xo), are the 
same as those of the case (b). 


The above two cases, (b) and (c), are grouped into the single case that <5 = 1/2 
and A‘^ = 2N, where is a positive integer. We depict the deflection angle for several 
values of N in Figure 2b. It should be noted that, unless N is unity, the deflection 
angle for small values of the impact parameter, h, can become a negative, increasing 
function of h in marked contrast with the case of the Schwarzschild spacetime. 

The appearance of a negative deflection angle indicates “reflection” of a light 
path, and we will later study under which conditions the light reflection occurs. 



12 


S. Tsuneishi, K. Watanabe and T. Tsuchida 


O' (deg) (a)N=1,2,3,4 «(deg) (b)N=1,2,3,4 




Fig. 2. The deflection angle, a{xo), is shown for some specific values of the parameters, {5, A), 
where is defined by Xo = l/cos^, and a{xo) is given in the unit of degree. The model 
parameters are (a) 5 = 0, = 2N and (b) S = 1/2, A'^ = 2N. We use a thin line for = 1, a 

short broken line for N = 2, a long broken line for N = 3 and a bold line for = 4. 


4.2. Optical scalars 


Since the local geometrical nature of the spacetime is described by the Riemann 
curvature, it is important to investigate the properties of optical scalars of the irro- 
tational null geodesic congruence.^ Let E( 2 ),E( 3 ),E( 4 )} = {k,rh,t,t} be a 

null tetrad such that 


= 0, + 5t6l (4-14) 


Then the optical scalars, namely, the expansion, 6, and the complex shear, u, are 
defined by 

9 = a = (4-15) 

The optical scalar equations are also conformally invariant, and we define the fol¬ 
lowing quantities associated with 


e = A^e-k^^V^\iiA, a = A^a, t>^ = Ai>^, = Ai>^. (4-16) 

Then the optical scalar equations becom^^’^ 
df) 1 

— + 0‘^ + |f7|^ = ——Rfj_yk^k^ = —TZ, (4T7) 

^ + 20a = = F, (4-18) 

where C^avg is the Weyl curvature. 

Dyer has obtained general solutions to the optical scalar equations in the general 
static, spherically symmetric spacetime.^ We summarize his results in Appendix C, 
and one finds that the Weyl source-term, E, and the shear, a, can be regarded as real 
quantities without loss of generality.^ Though the scalar-tensor-Weyl solution is 
not spherically symmetric, we can apply Dyer’s method as concerns the null geodesics 
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on the equatorial plane, y = oP By defining the new optical scalars, (7+ and C-, 
as 

— lnC'+ = 0 + (T, —lnC- = 9 — a, (4-19) 

dA dA 

Dyer has obtained the following equations, 


d^C+ 

dA 2 


i-n + F)C+, 


d^C_ 

dA 2 


{-n-F)C-. 


(4-20) 


The equations are equivalent to the geodesic deviation equations, and the 

variables, (7+ and (7_, represent, respectively, the size of the major axis and the 
minor axis of the infinitesimal image in the Einstein frame. We dehne the conformal 
transformation of C± by C± = AC± such that this transformation law is compatible 
with the others. 

The Ricci and Weyl source-terms, TZ{x) and E(x), in the Scalar-tensor-Weyl 
solution on the equatorial plane, y = 0, are calculated as (Appendix C, Ref.7) 


n{x) = 


Z \2 -^2 


X 


0-2 x2(x2 — 1) Yx2 — 1 

/\2 _ fjnr 




^ _ jh/crf f x-1 
x2 — 1 -I-1 


2(5' 


(x 2 — 1)2 

1 1 


(4-21) 


x3(x2 - 1} {x2 - 1 

F(x) = 6(A^ - 1) 

2(5 


X F{x), 


+ + (4-22) 

X J. \ X “1“ X / 


where h is the impact parameter as before. In general relativity, we have A = 5 and 
lZ{x) = 0. In the asymptotic region such that x » 1, lZ{x) and F{x), respectively, 
approximate to 


TZ{x) ~ 
F{x) ~ 


m2 (Z\2 - (52) 


^2^4 

3mh‘^ 


1 + 


4m /i2 


m\ 

1-1-j 

r / 


1 - 


A^-6^ 


m 


352 




4Z\2 -|- 52 — 5 777,2 
6(52 ^2 




Z \2 — 1 m^ 

(52 7 ^’ 


(4-23) 

(4-24) 


in terms of m/r, where x is related to r by with y = 0. 

It is immediately found that the Ricci source-term, TZ{x), satisfies TZ{x) > 0 and 
necessarily vanishes at the perihelion, x = Xq, because vanishes there. As for the 
Weyl source-term, F{x), its value is not always positive and can become negative 
along the light path. Since the Weyl source-term is mostly dominant over the Ricci 
source-term as will be found by the later numerical study, we classify the parameter 
space, {6, A), according to the qualitative properties of the Weyl source-term as 
follows. 
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When 5 = A = 1, the scalar-tensor-Weyl solution becomes the Schwarzschild 
spacetime, and we have 


“ih? 2,mh? 

(T^(x -|- 1)® r® 

One notes that the Weyl source-term has the following properties: 


(4-25) 


A1 It is strictly positive definite for x > Xq- 
A2 It is a monotonously decreasing function of x. 

A3 It takes the maximum value at the perihelion, x = xq, along the scattering orbit. 


Because of the properties above, one may simply expect that light rays are 
affected by the Weyl source-term mostly near the perihelion. For this reason, we 
first classify the parameter space, (d. A), according to the value of F{xo)- With (14-81) 
and ( 11 ^ . one finds that 


\ ^ 1 \ 

F(Xo) = — 3. 2 -TT ^- 7 

Xo'AXo^ - 1 ) \Xo^ - 1 / 

f{xo) = 36xo^ — + 26)xo + 2(5(A^ 


f{Xo), 

-!)• 


(4-26) 


The most intriguing thing is that F{xo) becomes zero for Xq = xf+ and Xq = xp-, 
where 

= ^ {(A^ + 26 ^) ± , 

D={A^ - - 24(5^ (2(5 - 1) {25 + 1). (4-27) 

Note that xf± and D are, respectively, roots and a discriminant of the quadratic 
equation, f{xo) = 0. Since Xq should be real and moreover satisfy the condition, 
Xo > 25 for 5 > 1/2 and Xq > 1 for (5 < 1/2, we classify the parameter space, {5, A), 
on the basis of the number of the conditional roots as follows: 


(a) The region of (5 > 1/2 

The discriminant, D, should be positive or zero in order for xf± to be real, and 
this condition is satisfied for the two cases, namely the case of A^ < A_^ and the 
case of A^ > A_|_^, where 

A| = 10(5^ ± V24(52 (2(5-l)(2(5-M). (4-28) 

For 5 > 1/2, A_|_^ and A_^ are always real and moreover satisfy the relation, 
0 < A _2 < A_|_^. In this region, a further condition, Xq > 25, should be satisfied 
for the scattering orbit, and we find that it holds only for the case, A^ > A+^, 
and is always violated for the case, A^ < A_^. We also hnd that the case of 
xf- < 25 < xf+ does not occur. The equations, A = A+ and 5 = 1/2, define the 
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critical lines on the parameter space, ((5, A), on which = xp- and Ap'^ = AJ^ 
hold, respectively. In summary, the region of (5 > 1/2 is classihed into the following 
two regions: 

(i) < Z\+2 

We shall refer to the region determined by the conditions, (5 > 1/2 and A?‘ < 
Z\+^, as the region Na (N: none) because the equation, F{xo) = 0, has no 
conditional root in this region. 

(ii) A’^ > Z\+^ 

We shall refer to the region determined by the conditions, 5 > 1/2 and A/^ > 
Z\+^, as the region II (II: two) because the equation, F{xo) = 0, has two condi¬ 
tional roots, XF+ and xp--, in this region. 


(b) The region of 0 < 5 < 1/2 

The quadratic equation, f{xo) = 0, always has two real roots, xpp and xp-, 
because its discriminant, H, is positive-definite for 0 < d < 1/2. In this region, a 
further condition, Xo > 1 , should be satisfied for the scattering orbit, and we find 
that xp- is inevitably less than unity. The condition, xf+ > 1, can be reduced 
to 5 < and we have xf+ = 1 on the new critical line defined by <5 = In 
summary, the region of 0 < d < 1/2 is classified into the following two regions: 

(i) 5 < 

We shall refer to the region determined by the conditions, 0 < <5 < 1/2 and 

5 < A^, as the region I (I: one) because the equation, F{xo) = 0, has a single 

conditional root, xpp, in this region. 

(ii) 5 > A^ 

We shall refer to the region determined by the conditions, 0 < d < 1/2 and 

5 > A^, as the region Nb because the equation, F{xo) = 0 , has no conditional 

root in this region. 


(c) The critical line, 5 = 1/2 

We have a single conditional root, xpp = 2(Z\^ — l)/3, for A^ > 5/2 and no root 
for < 5/2. 

(d) 5 = 0 

For any positive value of A‘^, F{xo) is negative. 


The results of the above classification are shown in Figure 3. 



F{x)/{2m) ^ F(x)/{2m)' 
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Fig. 3. We classify the parameter space, {S,A), into four regions (Na, Nb, I and II) according to 
the number of roots of a quadratic equation, f(xo) = 0, with the additional conditions described 
in the text. Horizontal and vertical axes denote 5 and A, respectively, and Z\ > 5 in the scalar- 
tensor-Weyl solution. Three critical lines, <5 = 1/2, S = A^ and A — Z\+, are shown, where A+ 
is given by A\ — 105^ + 5^24 (25 — 1) (25 + 1). Two specihc points, (5, A) = (0, 0) and (5, A) = 
(1,1), represent the Minkowski spacetime and the Schwarzschild spacetime, respectively. When 
we show some typical examples of the numerical results, we will mostly consider the points 
denoted by x as standard ones (see Section 5.1). 


(a) A=2.0, 6=0.7 (b) A=0.3, 6=0.3 



0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 1 1.2 1.4 1.6 1.8 2 


Fig. 4. The Weyl source-term, F{x), is shown as a function of x {xo < x) for several values of 
A and 5: (a) (5, Zl) = (0.7,2.0) in the region Na, (b) (S,A) = (0.3,0.3) in the region Nb, 
(c) (5, A) = (0.45, 0.8) in the region I and (d) (5, A) = (0.55, 2.5) in the region II. 
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In Figure 4, we show the Weyl source-term, F{x), as a function of x {xo < x) 
for several values of the parameters, 6 and A, and the impact parameter, h. While 
F{x) has similar properties as those in the case of the Schwarzschild spacetime for 
large h (see Al~A3), the following quite different properties B1~B3 can appear for 
small h. 

B1 It can be negative, and the maximum number of roots of the equation, F{x) = 0, 
is two for X > Xo- 

B2 It can have extremes, and the maximum number of the extremes is two for 

X > Xo- 

B3 It does not necessarily take the maximum value at the perihelion, x = Xo, along 
the scattering orbit. 

We numerically investigate the function, F{x), and classify the parameter space, 
{5, A), by the maximum number of the roots (Bl) and the maximum number of the 
extremes (B2). The results are summarized in Figure 5. In Figure 5a, the maximum 
number of the roots is shown by symbols, * and +. In the region denoted by 
there exists at least one impact parameter, h, such that F{x) = 0 has two roots. 
In the region denoted by -|-, F{x) = 0 cannot have two roots for any h, and there 
exists at least one impact parameter, h, such that F{x) = 0 has a single root. In 
the blank region, F{x) = 0 cannot have a root for any h- The maximum number of 
the extremes is shown in Figure 5b in the similar manner as in the case of Figure 
5a, namely * for two extremes, + for one extreme at the maximum and blank for no 
extreme. One immediately finds close resemblance between these classified regions 
in Figures 5a,b and the previous one in Figure 3. We will numerically show that 
the classification presented in Figure 3 really has a close relationship with several 
properties of light propagation and gravitational leasing. 


(a) (b) 




Fig. 5. We classify the parameter space, (<5, A), according to the qualitative properties of the Weyl 
source-term, Fix): (a) classification by the maximum number of roots of F{x) = 0 and (b) 
classification by the maximum number of extremes of Fix). Further details of these figures are 
explained in the text. 
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§5. Numerical results on null geodesics and optical scalars 


Now we will show several numerical results on null geodesics and optical scalars. 


5.1. Preliminaries 


There are three model parameters, a, 6 and A, in the scalar-tensor-Weyl solution, 
and we reduce the number of the parameters to two by hxing the values of m = 5a. 
In other words, m plays the role of a unit of length. When we show some typical 
examples of the numerical results, the following values of the model parameters, 
(5, Z\), will mostly be considered as standard ones (see Figure 3), 


{6, A) 


'[ha] (0.7,2.0) 
[hb] (0.3,0.3) 
[i] (0.43,2.0) 
Jii] (0.53,2.0) 


in the region Na 
in the region Nb 
in the region I 
in the region II 


(5-1) 


In Section 4.2 we have investigated the conditional roots of the equation, F{xo) = 
0, and find the single root, xf+, in the regions I and the two roots, xf+ and xf- in 
the region II. These roots, xf+ and xf-, are related to the impact parameters, hF+ 
and Hf-, by 

hF+ = I 

\XF+ - 1 

hF- = a^/xF^ - 1 , (5-2) 

\XF- - 1/ 



respectively. We will repeatedly emphasize the importance of these specific values 
of the impact parameter, hF+ and Hf-. 

On the equatorial plane, y = 0, a spatial part of the metric ra becomes 


2 / 1 \ —5 

, 9 m / X — 1 


- 1 


X^ 


A ^-1 


dx^ -|- {x^ — l)d(fP‘ 


(5-3) 


and the equatorial circumference, is given by 


Lx — 2fRx{x'), Rx{x') — ^ y _j_ 


\/x^ — 1. 


(5-4) 


On the basis of this geometrical consideration, we define the dimensionless, locally 
Cartesian-like coordinates, (X,Y), by 


2m 


2m 


(5-5) 
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for the sake of later convenience. 

The optical scalar equations should be solved along the corresponding light path, 
and therefore, we numerically integrate the geodesic equations and the optical scalar 
equations simultaneously under the initial conditions as follows: 


(a) The direction of the X-axis is chosen such that a light ray initially propagates 
in the opposite direction to the X-axis, namely from “right” to “left” in Figure 6. 


(b) The affine parameter. A, is chosen such that A = 0 at the initial time. 


(c) Let Xi be an initial value of the coordinate, x, of the null geodesic at A = 0. We 
should determine Xi such that an initial circumference radius, Ri = Rx{xi), is an 
approximately parameter-independent value. For Xi = 1005, we have 


Ri/m 


1 f xj-i y'^ 
5 Vxj + l) 


\lx^ -\ 


5 


i + - + o{ 

Xi 


i- 


101 O ( 10 " 2 ) 


(5-6) 


Unless otherwise noted, we let a value of Xi be 1005. 

Along the scattering orbit, x decreases from Xi to its minimum values, Xq, and 
then increases to Xi. When x becomes Xi again at A = A/, numerical integrations 
are terminated. We shall call A = A/ the final time. Under this geometrical 
configuration, we define the Einstein radius, /i^;, by 


hE 

2m 



(5-7) 


The Einstein radius, Le, usually defined in the weak gravity region such that 
m/hE « 1, and the above estimate of is compatible with the assumption that 
the gravity is comparatively weak even at the Einstein radius. 


(d) As for the optical scalers, their initial conditions at A = 0 are chosen as 


C+ = C-= 0, 


dC+ 

~d\ 


dC- 

~d\ 


such that we have (7+ = C- = A in the Minkowski spacetime. 


(5-8) 


Now we should remind ourselves that the optical scalars are conformally variant 
as C± = AC±. We shall therefore define a conformally invariant optical scalar, C, 
by 


_ C+ - C_ _ C+ - C_ 

~ c+ ~ c 


Since (7+ and C- represent the size of the major axis and the minor axis of the 
infinitesimal image, respectively, the new variable, C, is a measure of the image 
distortion rate. When either (7+ or (7_ vanishes at the final time for some value of 
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the impact parameter, h, we shall refer to this specific value of the impact parameter 
as a caustic point. It is important to note that we have two kinds of caustic points, 
namely the caustic point associated with C- and the one associated with (7+. At 
the caustic point, we have (7 ^ 1 for (7_ ^ 0 and |(7| ^ oo for (7+ ^ 0. 

5.2. Null geodesics and the deflection angle 

We show the trajectories of light rays on the (A, y)-plane in Figure 6. The model 
parameters, {5, A), are chosen to be the standard ones and the impact parameter, 
h, is chosen as (a) h/2m = 3.0 and (b) hl2m = 2.17, respectively. For the case 
(b), the light path in the case [nA] fails to be a scattering orbit, and therefore, we 
omit this case in Figure 6b. Figure 6c is an enlarged picture of the central part of 
Figure 6b. When the impact parameter, /i, is comparatively large as in the case of 
Figure 6a, the gravity acts attractively, and the effects of the scalar held and the 
non-sphericality remain quantitative. In contrast to that, Figure 6c shows us that 
the gravity can act repulsively for smaller values of h, as has already been predicted 
by the existence of the negative dehection angle in the previous section. 

The appearance of the repulsive trajectories explicitly proves that “rehection” 
of a light path can actually occurs, and we now investigate the conditions for 
that. By the simple geometrical consideration, one hnds that the situation of 


30 

20 

{a} h/2m=3.0 

' A=2.b,6=0.7 

A=2.0,8=0.53- 

A=2.0.8=0.43 - ■ 

A=0.3,6=0.3 - 

10 

" 

^ 0 


-10 

.X''' // 

■ / 

-20 

// 

-30 

, /. 


-40 -30 -20 -10 0 10 20 30 40 

X 

(b)h/2m=2.17 (c) h/2m=2.17 




Fig. 6. The trajectories of light rays are shown on the (X, y)-plane for (a) h/2m = 3.0 and (b) 
h/2m — 2.17. The enlarged central part of (b) is shown in (c). 










Light propagation and gravitational lensing on the Weyl-like spacetime 21 


4 

3.5 
3 

2.5 

< 2 

1.5 
1 

0.5 
0 

0 0.5 1 1.5 2 

5 

Fig. 7. The reflecting region is shown by symbols, +, on the classified parameter space, {S,A). 



d/dX{dY/dX) > 0 inevitably occurs along a reflecting, counterclockwise light path 
around the singularity at x = 1. We succeed in reducing this condition to the 
occurrence of p{x) < 0 , where an auxiliary function, p(x), is defined by 


PM = ^4^3^71)2 (^) XlW- 

q(x) = {{S^ - A^)x + {A‘‘ - 1)(5} (Itt) “ 

+'■'"(1^1) 



(5-10) 


In particular, we have p{x) = “imh?/r^ > 0 in the Schwarzschild spacetime. 

In Figure 7, we numerically compare the previously classified four regions, Na ~ 
II, with the reflecting region in which the inequality, p{x) < 0, can apply. Though 
their relationship is halfway close, it is an intriguing thing that some parts of the 
boundary of the reflecting region coincide with the critical line, 5 = 1 / 2 . 

As for the deflection angle, we numerically integrate with (EIHI) and depict 
it as a function of the impact parameter, h, in Figure 8 , where the marks, /ii, /i 2 and 
/ 13 , represent the following specific values of the impact parameters, respectively: 

(/ii, / 12 ) = {hp-, hF+) in the case [ h ], ^3 = ^f+ in the case [ i ]. (5T1) 

We numerically find that the deflection angle generically has two extremes near 
h = hpp and h = hp- in the region II and a single extreme near h = hpp in the 
region I. Graphs of the deflection angle in the regions, I and Nb, have the similar 
properties as follows: 


Cl It has a single extreme and the extreme is the maximum. 

C2 The deflection angle remains to be finite as the impact parameter approaches 
to zero. 
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h/2m 

Fig. 8. The deflection angle, a, is shown as a function of the impact parameter, h. The model 
parameters are chosen to be the standard ones. 


(a) h/2m=3.0 



(b) h/2m=2.17 



Fig. 9. The deflection angle, a, is shown as a function of the parameter, A, for (a) h/2m = 3.0 
and (b) h/2m = 2.17. Thin lines show the deflection angles in the Weyl solution (A = 5). For 
general cases of > 5, lines corresponding to the constant 5 have terminal points (x in the 
fignres) on the thin line. Open circles and filled circles denote the data points snch that Q = 1 
and Q = —1, respectively, where Q = {Ar — l)/5^ is a measure of the non-sphericality of the 
spacetime symmetry. An open triangle in (b) denotes the case that the scattering orbit cannot 
exist. 


In particular, the properties, Cl, are unexpected ones because of the absence 
of the characteristic scales, and hp-, in the region Nb- Let ha be a new 

characteristic value of h at which a graph of the deflection angle has the maximum. 
We numerical find that ha continuously changes beyond the boundaries between the 
region I and the region Nb and that ha = 2m ~ 6m in the region Nb- We also find 
that ha is a decreasing function of S and is an increasing function of A. 

By fixing the values of 6 and h, we show the deflection angle as a function of A 
in Figure 9 and find that the deflection angle is always a decreasing function of A. 
We also show the cases in general relativity {A = S) for comparison. 

By further numerical investigation on the parameter space, (5, Ll), we find the 
following notable properties of the deflection angle, a{h). 

D1 While it is positive definite in the Schwarzschild spacetime, it can be negative 
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in the Scalar-tensor-Weyl solution. 

D2 While it is a monotonously decreasing function of h in the Schwarzschild space- 
time, it can have two extremes at the maximum. 

D3 Since a circular orbit exists when 5 > 1/2, a light ray passing by the circular 
orbit can wind around it many times, and the deflection angle becomes very large. 
In contrast to the cases of 5 > 1/2, it always remains to be finite in the case of 0 < 
5 <1/2 even when h approaches zero. In the case of (5 = 1/2, we analytically find 
that it remains to be finite for Z\ > 1 and diverges for A <1 when h approaches 
4m. 

We classify the parameter space, (5, A), according to the number of extremes of 
a{h) and compare this new classification with the previous one, as is shown in Figure 

10. The similar classification according to the occurrence of negative deflection angles 
is also done, as is shown in Figure 11. One may immediately note the surprising 
coincidences among these results found in Figures 10 and 11. Though we have not 
yet succeeded in explaining the reasons for these coincidences, we expect that the 
classification shown in Figure 3 may have some physical significance for gravitational 
lensing properties. As for the light reflection and the deflection angle, we find that 
their properties in the region Nb is quite similar to those in the region I. 

Finally, it should be noted that this kind of the light reflection can be possible 
even in general relativity. That is, the light reflection inevitably occurs even on the 
line, A = 5, for 5 < 1/2. On the other hand, the appearance of the region II must 
be an intrinsic effect to the scalar-tensor theories of gravity, and we expect that the 
Il-like region may always exist not only in the scalar-tensor-Weyl solution but also 
in most of the other solutions. While the boundary between the NA-like region and 
the I-like region will always be clear, the boundary between the NA-like region and 
the Il-like region may, however, be slightly vague, as has seen in Figures 7, 10 and 

11 . 



Fig. 10. We classify the parameter space, 
((5, Z\), into three regions according to the 
nnmber of extremes of a(h): * for two ex¬ 
tremes, -I- for a single extreme and blank for 
no extreme. 



Fig. 11. We classify the parameter space, 
((5, Z\), into two regions according to the oc¬ 
currence of negative deflection angles. In the 
region with filled circles, a deflection angle 
can be negative. 
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5.3. Optical scalars and the image distortion rate 

We numerically solve the optical scalar equations along the light path by the 
procedure described in Section 5.1. In Figure 12, we show (a) C+ and (b) C- as 
functions of the affine parameter, A, for several values of the impact parameters, h. 
The model parameters are chosen as those of the case [nA], where ((5, A) = (0.7, 2.0). 
One must find at first sight that gravitational effects of the source-terms, F and IZ, 
on the light rays appear as a rapidly changing slope of the graph of (7+ in the vicinity 
of A = 50. As for C-, the slope changes once for /i/2m =3.0 and 4.0 and twice for 
h/2m = 2.5. This “bending” of the graph is found in the narrow range of A, and 
its width, Z\A, is 2 ~ 3 for the case of C- with hj^m = 2.5 and invisibly thin for 
the other ones. In this sense, a picture of the thin lens is applicable, and, roughly 
speaking, the gravitational effects appear as a bending angle of the graph. 


(a) A=2.0,5=0.7 (b) A=2.0,5=0.7 




Fig. 12. We show (a) C+ and (b) C- as functions of A for several values of the impact parameters, 
h, as denoted in each figure. The model parameters are chosen as those of the case [nA], where 
(5,A) = (0.7,2.0). 

In Figures 13, 14 and 15, we show (a) (7+ and (c) (7_ for the standard cases, 
[ub], [ i ] and [ h ]. For comparison, the enlarged pictures of these graphs are also 
shown in (b) and (d) for (7+ and (7_, respectively. In any case, a rough picture of 
the thin lens is applicable again. To confirm the forecast that the lenses are thin, 
we show the Ricci and Weyl source-terms as functions of the affine parameter. A, in 
Figures 16~19. These thin lenses are, however, divided into the following two groups 
by their thinness (or thickness). 

El A lens is invisibly thin such that the bending angle of the graph is a moderately 
varying function of the impact parameter, h. 

E2 A lens has a visibly finite thickness, AA, such that one recognizes a comparatively 
rapid variation of the graph when light rays goes through the lens. 

It should be noted that the above classification of lenses is brought not so much 
by the intrinsic properties of each classified region in the parameter space as by the 
extrinsic factors, especially the impact parameter, h, and the initial distance, Ri, as 
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is explained by the following simple model. 

Under the thin lens assumption, we approximate the optical scalar equations 
(I4-2UI) by 

d'^u± _ [o (|A|>e) 

dX'^ (|A| < e), 

where and e ~ AXj2 are an amplitude of the source-term and a thickness of the 
lens, respectively. The affine parameter. A, is redehned such that the center of the 
thin lens is at A = 0. For example, when the combined source-term, —TZ + F, of (7+ 
is negative, the negative sign of (I5T2I) is chosen for C^: 


d'^C. 


+ 


= (-7^ + F) C 4 


d^u- _ 10 (|A| > e) 

I —K^u- (|A| < e). 


dA2 ' ' ' ^ dA2 

The equation is easily solved under the initial conditions, 

du± 


u± = 0 , 


dX 


= 1 at A = —Ao < 0, 


(5-13) 


(5-14) 


and we find the following properties under the condition that e ~ AXj2 << Aq ~ 


FI The case of CK << 1 

It corresponds to the case of El. The bending angle, A9±, of the graph of C± is 
well approximated as 

Ae± = ±{kXo){2k€) ~ Fk'^XoAX. (5-15) 

Note that \A0± \ can be larger than unity for the cases of large values of kAq ~ nRi. 
F2 The case oi en > 1 

It corresponds to the case of E2. When light rays pass through the lens,ri_ oscil¬ 
lates approximately eft/vr times, while exponentially grows as 



log u+ kX + log 


(5-16) 






C-/2m C+/2m CW C-/2m C+/2m 
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(a) A=0.3,5=0.3 


(b) A=0.3,6=0.3 




W2m W2m 


(c) A=0.3,6=0.3 


(d) A=0.3,5=0.3 




13. We show (a) C+ and (c) C_ as functions of A for several values of the impact parameters, 
h, as denoted in each figure. The model parameters are chosen as those of the case [ub], where 
{S,A) = (0.3, 0.3). The figures (b) and (d) are the enlarged parts of the figures (a) and (c), 
respectively. 


(a) A=2.0,5=0.43 (b) A=2.0,6=0.43 




(c) A=2.0,5=0.43 (d) A=2.0,6=0.43 




Fig. 14. We show (a) C+ and (c) C- as functions of A for several values of the impact parameters, 
h, as denoted in each figure. The model parameters are chosen as those of the case [ i ], where 
{S,A) — (0.43,2.0). The figures (b) and (d) are the enlarged parts of the figures (a) and (c), 
respectively. 
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(a) A=2.0,5=0.53 (b) A=2.0,8=0.53 




(c) A=2.0,5=0.53 


(d) A=2.0,5=0.53 




W2m 


X/2m 


Fig. 15. We show (a) C+ and (c) C- as functions of A for several values of the impact parameters, 
h, as denoted in each figure. The model parameters are chosen as those of the case [ ii ], where 
{S,A) — (0.53,2.0). The figures (b) and (d) are the enlarged parts of the figures (a) and (c), 
respectively. 
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(a) A=2.0, 5=0.7 (b) A=2.0, 5=0.7 




(c) A=2.0, 5=0.7 


(d) A=2.0, 5=0.7 




16. We show the Weyl and Ricci source-terms, respectively, in the figures (a) and (c), as 
functions of the affine parameter, A, for several values of the impact parameter, h. The model 
parameters are chosen as those of the case [nA], where ((^, ^) = (0.7, 2.0). The figures (b) and 
(d) are the enlarged central parts of the figures (a) and (c), respectively. 


(a) A=0.3, 5=0.3 


(c) A=0.3, 5=0.3 




(b) A=0.3, 5=0.3 



5 h 
4 

3 h 


(d) A=0.3, 5=0.3 


h/2m=3.0 • 
h/2m=2.17 • 
h/2m=1.0 • 
h/2m=0.5 • 


51 

W2m 


Fig. 17. We show the Weyl and Ricci source-terms, respectively, in the figures (a) and (c), as 
functions of the affine parameter. A, for several values of the impact parameter, h. The model 
parameters are chosen as those of the case [ne]: = (0.3, 0.3). The figures (b) and (d) are 

the enlarged parts of the figures (a) and (c), respectively. The Ricci source-term, IZ, is exactly 
zero for any value of h. 
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(b) A=2.0, 5=0.43 
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18. We show the Weyl and Ricci source-terms, respectively, in the figures (a) and (c), as 
functions of the affine parameter, A, for several values of the impact parameter, h. The model 
parameters are chosen as those of the case [ i ], where (<5, Z\) = (2.0,0.43). The figures (b) and 
(d) are the enlarged parts of the figures (a) and (c), respectively. 


(a) A=2.0, 5=0.53 (b) A=2.0, 5=0.53 




(c) A=2.0, 5=0.53 
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(d) A=2.0, 5=0.53 
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Fig. 19. We show the Weyl and Ricci source-terms, respectively, in the figures (a) and (c), as 
functions of the affine parameter, A, for several values of the impact parameter, h. The model 
parameters are chosen as those of the case [ h ], where (J, A) = (0.53, 2.0). The figures (b) and 
(d) are the enlarged parts of the figures (a) and (c), respectively. 
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The oscillatory properties of C- have been expected by Dyer’s results (see Ap¬ 
pendix C), while we find that C+ can oscillate because its combined source-term, 
—TZ + F, can be negative not only in the scalar-tensor-Weyl solution but even in 
Voorhees’s one. We have tried to find the numerical examples in which (7+ and 
C- oscillate. After finely adjusting the model parameters, (5,/\), and the impact 
parameter, h, we find several examples in which C- oscillates. As for (7+, however, 
we have failed. In Figure 20, we show the optical scalars, C±, together with their 
source-terms for the case [ ii ] with h/2m = 2.1609. With Figures 20a and 20b, we 
obtain e = 0.04 x (2m) and k = 77/(2m) such that en = 3.1 > 1. We then estimate 
the expected number of times of the oscillations at ck/tt ~ 1, which is in accordance 
with the numerical one (see Figure 20f). As for (7+, the analytic result (I5T6I) is 
approximate to the numerical one as follows. 


log (7+/2m 


78A/(2m) -F4 
80A/(2m) -h 6 


(Analytic) 

(Numerical). 


(5-17) 


When the combined source-term has plural peaks, one can superpose the correspond¬ 
ing thin lenses. 


Now we summarize the qualitative properties of C+ and (7_ along the light path 

for each region. 

(Na) In the region Na, the Weyl source-term, T, is always dominant over the Ricci 
source-term, IZ, and moreover is positive-definite. When the impact parameter is 
larger enough than m, the inequality, en « 1, holds, and the bending angle of the 
graph monotonously increases as the impact parameter decreases. Accordingly, 
there is a single caustic point, h = ha, near the Einstein radius, He, such that 
C- vanishes at the final time, X = Xf. For the smaller impact parameter, the 
inequality, eK > 1, holds, and (7_ can oscillate. 

(Nb) In the region Nb, the Weyl source-term, F, is always dominant over the Ricci 
source-term, TZ. When the impact parameter is large enough, the properties of (7+ 
and C- are similar to those in the region Na. For the smaller impact parameter, 
the Weyl source-term can be negative such that (7+ can decrease, while (7_ can 
increase. Accordingly, there is a new caustic point, h = hf^, with = 3m ~ 6m 
such that (7+ vanishes at the final time. The intriguing thing for h < hf, is that (7+ 
has negative values with the very large amplitude, while (7_ is moderately positive. 
The case of (7_ with hj2m = 0.5 may be an example of the superposition of triple 
lenses (see Figure 17a). 

( I ) In the region I, the Weyl source-term, F, is mostly dominant over the Ricci 
source-term, TZ, and the Ricci-dominant cases are found only in the very narrow 
range of the impact parameter around hF+- When hF+ < hE, a first caustic point 
exists near He, and a second caustic points, h = h},, appears near hF+ such that 
(7+ vanishes at the final time. It should be noted that the inequality, hF+ < hE, 
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(a) h/2m=2.1609 


(b) h/2m=2.1609 
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(f) h/2m=2.1609 



Fig. 20. In the case [ ii ], where (<5, A) = (0.53, 2.0), we show (a) —TZ + F, (b) —TZ — F, (c) C+, 
(d) C- and (e) log(7+. The impact parameter, h, is finely adjusted as h/2m = 2.1609 such that 
C- oscillates. The figure (f) is the enlarged part of the hgure (d). 


is easily violated in our numerical studies with Xi = 100(5 (see Figure 23). In these 
cases, the caustic point, h = hb, becomes the hrst one. Let 6*+ and C- at the final 
time, A = Aj-, be functions of the impact parameter, h. For h < hb, C+ becomes 
negative with a larger amplitude for smaller h, while C- changes its sign again 
and becomes positively larger for smaller h. Accordingly, a third caustic point can 
exist such that C- vanishes at the final time. 

(II) In the region II, the Weyl source-term, F, is mostly dominant over the Ricci 
source-term, TZ, and the Ricci-dominant cases are found only in the very narrow 
ranges of the impact parameter around and hp-. We find the very compli¬ 
cated structure of caustic points. For example, in the case [ ii ], we find four caustic 
points, ha, hb, he and hd, where ha > hb > he > hd holds. The first caustic point, 
ha, appears near the Einstein radius, hp, such that C- vanishes at the final time. 
The second and third ones, hb and he, appear near hp+ and hp-, respectively. 
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Fig. 21. In the figure (a), the image distortion rate, C, is shown as a function of the impact 
parameter, h. The model parameters are chosen to be the standard ones. The figure (b) is the 
enlarged part of the figure (a). 


such that (7+ vanishes at the final time. Finally, the fourth one is associated with 
C- such that C- vanishes at the final time. 

We have defined a conformally invariant variable, C, by C* = ((7+ — (7_)/(7+, 
which can be interpreted as the image distortion rate. In Figure 21, we show C at 
the final time as a function of the impact parameter, h. As is noted above, one of the 
most intriguing properties of the graph is the appearance of the caustic points, and 
we find the following four cases of the number, A^+, of the caustic points associated 
with (7+: 

(a) N+ = 0, (b) iV+ = 1, (c) N+ = 2, (d) N+ / 0 and He < hF+- 

Now we should remind ourselves that He is an extrinsic scale including the initial 
distance to the source, Ri = Rx{xi). It is naturally expected that the cases of (d) 
are distributed to (b) and (c). As such an example, we investigate the cases of 
(5, A) = (2.0,0.3) in the region I and (5, A) = (4.0,0.6) in the region II. In both 
cases, we have hF+/‘2^m > hE/2m ~ 7 for Xi = lOOd and hF+f^m < hE/^m ~ 14 
for Xi = 4005. We show the image distortion rate, C, as a function of the impact 
parameter, h, in Figure 22. The graphs show the “anomalies” for the cases of 
Xi = 100(5 in the sense that the first caustic point, (7 = 1, near the Einstein radius, 
/iB, disappears. That is, a naked caustic point associated with (7+ appears near 
hF+- For the cases of Xi = 400(5, this anomalous caustic point is hidden by the 
Einstein radius, He, such that the first caustic point appears near He sts the one 
associated with (7_. We classify the parameter space, ((5, A), on the basis of the above 
conditions, (a)~(d), as is depicted in Figure 23. Again, we find a close relationship 
between the above new classification and the previous one. 

We summarize the relationship between our numerical results and the classifi¬ 
cation of the parameter space, ((5, A), in Table I. 















Light propagation and gravitational tensing on the Weyl-like spacetime 33 


(a) A=4.0,5=0.6 


{b} A=4.0,5=0.6 
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(c) A=2.0,5=0.3 


(d) A=2.0,5=0.3 




Fig. 22. We show the image distortion rate, C, as a function of the impact parameter, h, for two 
cases: (a) (5, A) = (0.6, 4.0) in the region II and (c) (5, A) = (0.3, 2.0) in the region I. The 
enlarged figures (b) and (d) correspond to the fignres (a) and (c), respectively. 



Fig. 23. We classify the parameter space, (<5, A), into the fonr groups, (a)~(d), on the basis of the 
number of the caustic points, A^+, associated with C+, as is explained in the text: (a) an open 
circle, (b) +, (c) and (d) a reversed triangle. 


Table I. The classification of the parameter space, {S, A) 


Classihcation of F{xo) 

Na 

Ns 

I 

II 

Fig. 3 

Reflecting region 

mostly blank 

+ 

+ 

+ 

Fig. 7 

Extremes of the deflection angle 

mostly blank 

+ 

+ 

* 

Fig. 10 

Negative deflection angle 

blank 

• 

• 

mostly • 

Fig. 11 

Caustic points associated with C+ 

mostly o 

+ 

+ 

* 

Fig. 23 
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§6. A closing remark 


In this paper, we have investigated several properties of null geodesics and null 
geodesic congruences with the aim of finding the crucial, experimentally detectable 
difference between the theory of general relativity and scalar-tensor theories of grav¬ 
ity. Our studies are very restricted ones because they are performed only for specific 
null geodesics (ones on the equatorial plane, y = 0) with a specific solution (the 
scalar-tensor-Weyl solution). These specihc null geodesics may be unstable when 
there are additional small fluctuations, especially those in the y-direction. One may 
doubt the realities of the specific solution as a mathematical model of the spacetime. 
Indeed we cannot guarantee its realities. Nevertheless, we would like to claim that 
we have succeeded in discovering the new world, namely the region II in the clas¬ 
sified parameter space. The leasing properties found in this region are never seen 
in the other regions, Na, Nb and I. Unfortunately, the border between the region 
II (the new world) and the region Na (the Einstein-like world) is rather vague, and 
moreover, the most precious things and hp-) may be hidden by the Einstein 

radius, hp. But one will go there: someday, somewhere, somehow^^ 

Acknowledgements We would like to thank Ken-ichi Oohara and Ryoichi Nishi 
for helpful comments. 


Appendix A 

The Ernst-like equations in scalar-tensor theories 


In the stationary, axisymmetric spacetime, there exist two Killing vectors, ^ = dt 
and rj = dff). It is shown that the energy-momentum tensor, of the Maxwell 

fields satisfies the following relations,^ 

= 0. (A-l) 

One finds that the similar relations hold for = ‘^d^ipdy^p, accordingly, with 
Theorem 7.1.1 in Ref. 13, the metric of the Einstein-Maxwell system with a scalar 
field, if, can be reduced to the following form, 

= —e^'^{dt — u>d4>)'^ e~‘^^ [e^'^{dp‘^ -\- dz^) -\- p^d(j)^] , (A-2) 

where the metric functions, V'j 7 and uj, are functions oi = p and = z. 

After long and complicated calculations, one finds that the field equations are 
reduced to the following set of equations, 

e^’/' = ^ (T + T) + m, 

= (VS) ■ [(VT) + 23V<P] , 

= (Vd>) • [(VT) + 2^V<?] , 


(A-3a) 

(A-3b) 

(A-3c) 
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7.1 = P [(V’,!)^ - (^^, 2 )^] - ^ [(w 1 )^ - (w 2 )^] 

- ^,2^,2)e-"’^ + p , (A-3d) 

7.2 = 2pV'7V',2 - - />(^ 7^,2 + ^, 2 ^,i)e“^^ + 2p(p^iip^2, (A-3e) 

Ip 

VV = 0 , (A-3f) 


where, for any functions, / and h, 


vV 


I d f df\ 52 / dfdh dfdh 


/,i 


' 


(A-4) 


In Ref.l4, there are explicit forms of the metric function, cu, and the Maxwell 
fields in terms of the complex potentials, £ and L>. One notes that the equations 
(A.3a)~(A.3c) are the same as those in general relativity and are referred to as the 
Ernst equations of the Einstein-Maxwell system. The remaining equations contain 
the scalar field contributions, and, as for the metric functions, effects of the scalar 
field appear only in 7 . 


Appendix B 

The scalar-tensor-Weyl solution 


In the case of static, axisymmetric vacuum solutions, one has oj = = 0, and 

the field equations (IA-3al) (Toa become the following set of equations, 

e^^ = £, (B-la) 

£\/^£ = (VT) • (VT), (B-lb) 

7.1 = P - (V’, 2 )^] + p - iT,2?] , (BTc) 

7.2 = 2pV',iV’,2 + 2pip^np^2, (B-ld) 

= 0, (BTe) 

where is a real function. One finds that the equation for '0 is reduced to 

= 0. (B-2) 

That is, ip and 1 /^ are harmonic functions. 

We introduce oblate and prolate coordinates, {x,y), defined by 

P = fT^(a:2 + e)(l-i/2)^ (B-3a) 

z = crxy, (B'3b) 

where cr is a positive constant, and e = ±1. The cases, e = 1 and e = —1, are referred 


to as oblate and prolate, respectively. The equation (IFl) is then reduced to 
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A similar equation holds for ip. 

The simplest prolate solution for ^ is given by 


Ip = — In 


X — 1 


^x -h 1 / ’ 

where S is an integration constant. A similar solution for (p is given by 


d /x — 1 
ip = <Po + -In 


X -h 1 / ’ 


(B-5) 


(B- 6 ) 


where ipo and d are integration constants. Then the corresponding metric function, 
7 , becomes 


e2^ = 


2 1 N 

X — 1 




where 

Z\2 = 6^ + d^. 

For completeness, we give an explicit form of the metric: 


(B-7) 

(B- 8 ) 


ds^ = — 


X — 1 

X + 1 


dt^ + cj^ 


X — 1 


-5 


x^ — 1 

x2 _ y2 




(x^ - y^) 


X + 1^ 

dx^ dy"^ 


+ 


+ - 1)(1 - V^)d'f 


(B-9) 


x^ — 1 1 — y^ 

One notes that, though xp and ip are independent of y, 7 depends on both x and y. 

Appendix C 

- Analytic solutions to the optical scalar equations - 


C.l. Static, spherically symmetric spacetime 


We summarize the analytic results obtained by Dyer.^ A metric of the static, 
spherically symmetric spacetime is given by 

ds^ = -e^^df + e^^dr^ + e‘^^df2^, (CT) 

where A, B and C are functions of r. A null vector, /c^, which is tangent to the null 
geodesic in this spacetime, is obtained as 

fe0 = e-2^, A:2 = 0, k^ = he-^^, 

’_ (C-2) 

= ±e-(^+^Vl-h2e2(^-^), 

where a constant, h, is an impact parameter, and we assume that the geodesic is on 
the equatorial plane, 9 = 7 r/ 2 , without loss of generality. The complex null vector. 
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t^, becomes 


t^ 


S± 


^-c 


2 V 2 [h 


S+S-e^ + H{Si + Si) 




-{Si + Si)e^ + 2HS+S- 


2 V 2 [h 
■ 'J e“*" ± he~^, H = —^ J 


‘ e-^ 


B' 


S+S- 


i 

Hh 

V 2 ' 

^-^dr, 


t^ = 


-2B 


(C-3) 


where a prime denotes a differentiation with respect to r. Then the Ricci and Weyl 
source-terms, IZ and F, in 6311 and are evaluated as 

-TZ + F = e-2(^+^) [B" + {B')^ - B'C' - B'A') 

- /i2e-2(^+-B) (^c" + (C')^ - C'B' - C'A') , 

-n-F = e-2(^+^) [B" + (R')^ - B'C' - B'A') 

- /i2e“2(^+^) ( qh _ _ 


Since the Weyl source-term, F, is real, one can let the shear, a, be real without loss 
of generality. Dyer has introduced the new optical scalars, C±, defined by 

^ In C± = 0±cT, (C-5) 

dX 

and has obtained the following equations, 

^ = (-F±F)C±. (C- 6 ) 

General solutions to (|G- 6 ll are given as 


C+ = G°\/e2-B - h'^i 


C- = C^e^sin ( h 


.2C 


{/' 


^A+B+C 


^q2B _ /j2g2Cj 2 


-dr + D 


+ (’ 


^q2B _ /j2g2C 


dr + Di 


(C-7) 


where and D% are integration constants. 


C. 2 . Static, axisymmetric spacetime 

A metric of the static, axisymmetric spacetime is given by 

ds^ = -e^^dt^ + e^^dx"^ + e^^ V + (C- 8 ) 

where the metric functions, a, (3, 7 and pL are functions of x and y. A null vector, 
k^, which is tangent to the null geodesic on the equatorial plane, y = 0 , is obtained 


as 


fcO = e = 0 , k^ = he 

k^ = ±e“^\/e“2“ — h‘^e~‘^f^, 


(C-9) 
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where a constant, h, is an impact parameter. The complex null vector, becomes 




1 


2^/2 

1 

2^/2 


^e^5+5_ + H{Sl + Sl) 

hsl + + 2HS+S. 

h 





S± = \/e^^^^^”/ie^. 



j-x a 

J 5 + 5 _^ 



—=Hhe 

V2 


2fi 

5 


(CTO) 
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